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We report measurements of diagonal susceptibilities for the baryon number, xb, electrical charge, 
Xq, third component of isospin, xiy strangeness, xs, and hypercharge, xy, as well as the off-diagonal 
XsQi Xby , Xbs, etc.. We show that the ratios of susceptibilities in the high temperature phase 
are robust variables, independent of lattice spacing, and therefore give predictions for experiments. 
We also investigate strangeness production and flavour symmetry breaking matrix elements at finite 
temperature. Finally, we present evidence that in the high temperature phase of QCD the different 
flavour quantum numbers are excited in linkages which are exactly the same as one expects from 
l/^ ■ quarks. We present some investigations of these quark-like quasi particles. 
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I. INTRODUCTION 



■ Experiments plan to study fluctuations of conserved quantities in heavy-ion collisions at the RHIC and LHC in 
^ ' different rapidity windows. With proper particle identification, one can measure in the experiment both absolutely 

, conserved quantities like the baryon number (B) and the electrical charge (Q), as well as quantities which are 

■ conserved only under the strong interactions, such as the third component of isospin (/a), the strangeness (5) and 
the hypercharge (Y). These observations can be used to extract fluctuations in the numbers of these quantities 
Such observations need to be compared to predictions of quark number susceptibilities (QNS) from lattice QCD. In 
this paper we report on lattice computations of a variety of diagonal QNS — xb, XQ: Xi^ Xs and xy- One of the main 

f*^ • results in this paper is the extraction of predictions for the ratios of these susceptibilities which survive the continuum 
^ ^ limit. Our second important result is the investigation of the strange quark sector of the theory: we extract the 
. Wroblewski parameter in a dynamical QCD computation for the first time, and also investigate the dynamics and 
I ■ kinematics of flavour symmetry breaking in QCD. Further, we present results on the cross correlations xbq, Xby, Xbs 
Oh, and XQY- These cross correlations are used to explore the charge and baryon number of objects that carry flavour. 
^ ' We find that the baryon number of flavour carrying objects immediately above the QCD crossover temperature, Tc, 
are 1/3 and the charges are 1/3 or 2/3. We find furthermore, that these objects are almost pure flavour — anything 
carrying u flavour has only tiny admixtures of d and s flavours, etc.. This is our third main result. 

We have bypassed the necessity of numerically taking the continuum limit of the theory by restricting attention 
to the high temperature phase where it is easy to define robust observables which have little, or no, lattice spacing 
' dependence. We demonstrate the robustness of the observables in quenched QCD, and then compute these quantities 
in QCD with two flavours of light dynamical quarks. These are also good observables in the sense of — 

XL <ji 

where xk and xl are QNS for the conserved quantum numbers K and L and aK and are the variances. The two 
variances must be obtained under identical experimental conditions, after removing counting (Poisson) fluctuations 
as suggested by 3]. Thus the robust lattice observables give predictions for robust experimental observables. 

Either oi K or L can also stand for a composite label (M, N) where M and N are conserved quantum numbers — 
in this case the susceptibility is an off-diagonal susceptibility, and the variance has to be replaced by the covariance 
of M and N. Note the relation with the correlation coefficient — 
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where again, the expressions are robust both on the lattice and in experiment. The study of these robust variables 
tells us about the relative magnitudes of fluctuations in different quantum numbers. The study of these quantities is 
one of the main results reported here. 

We further present investigations of the strange quark sector of the theory. The robust variable Csu is closely 
related to the Wroblewski parameter which can be extracted from experiments. This shows strong dependence on 
the actual strange quark mass, TOs, in the vicinity of Tc- Since rus ~ Tc, it seems that part of this sensitivity could 
be attributed purely to kinematics. We investigate the dynamical matrix elements which are responsible for flavour 
symmetry breaking in QCD and compare the importance of kinematics and dynamics in the strange quark sector. 
This is our second major result. 

One outstanding question about the high temperature phase of QCD is the nature of flavoured excitations. There is 
ample evidence that quarks are liberated at sufficiently high temperature — the continuum limit of lattice computations 
of screening masses are consistent with the existence of such a Fermi gas for T > 2Tc (Ji |5| ; quantitative agreement 
between weak coupling estimates of the susceptibilities 0,13 and the lattice data~0,|3| also confirm this; the equation 
of state at very high temperature also testifies to this. However, comparison of lattice results and weak coupling 
computations of these quantities fail for T < 2Tc. Our third new result concerns this matter of the thermodynamically 
important single particle excitations. 

We address this question in the most direct way possible — create an excitation with one quantum number and 
observe what other quantum numbers it carries. Technically, this involves the measurement of robust ratios of off- 
diagonal QNS; the correlation between quantum numbers K and L can be studied through the ratio 

_ {KL)-{K){L) 

C^KL)/L - ^i2)_(^)2 • (3) 

We find that such measurements are feasible on the lattice, and are open to direct interpretation. We also suggest 
that they could be performed in heavy-ion experiments, as direct tests of whether quarks exist in the hot and dense 
matter inside the fireball. A recent suggestion of ^] is the measurement of just such a variable: essentially C(^bs)/S- 

We find that, immediately above Tc, the baryon number, charge and other flavour quantum numbers are linked 
with each other in exactly the same way as they are in quarks. For example, excitations which carry unit strangeness 
carry baryon number of —1/3 and charge of -1-1/3. This, together with the fact that there is also a failure of weak 
coupling theory, would imply that the QCD plasma phase is a "quark liquid" in the sense that the quasi-particles 
carry the quantum numbers of quarks, but the interactions between them are too strong for the system to be treated 
in weak coupling theory. Extension of these measurements to finite chemical potential for T > Tc and fi ^ T could 
allow us to check whether or not the system is a normal Fermi liquid ITdi. Such an extension is feasible since the 
Taylor series expansion of the free energy in /i/T has a radius of convergence much higher than unity for T > Tc- 

This is an appropriate place to remark upon a few aspects of our computations. Having removed most of the lattice 
spacing uncertainties by using robust variables, we have to control only the quark masses. We do this partly by 
performing the computations in an approximation called partial quenching. In this approximation the valence quark 
masses in the theory are tuned keeping the sea quark masses fixed. We explore the dependence of the robust variables 
on the sea quark masses and find that the results are not very sensitive to these parameters. This is expected — away 
from a phase transition there is no more than a 5% change in the QNS in going from quenched to Nf = 2 dynamical 
QCD, and one expects the change to be smaller in going from Nf = 2 to Nf = 2 -I- 1, as long as one avoids the vicinity 
of the phase transition. The ratios are even less sensitive to the sea quark content than the QNS. In this study we 
have concentrated on the numerically more important effect of the valence quark masses. 

We have used two flavours of dynamical sea quarks of bare mass m = 0.1 to study a temperature range upto 
about 2Tc. These quark masses are such that nip/Tc = 5.4 and m-^/mp = 0.3 — which makes this the smallest quark 
mass used in a systematic study of fluctuations. We have taken the strange quark to be quenched and to have a 
bare mass in the range m^/Tc = 0.75-1. This gives the correct physical values of the ratio m^c/mp. We have also 
investigated the effect of decreasing the valence light quark mass by a factor of three in order to get at the same time 
the correct physical value of the ratio rriT^/mp, and varying the strange quark mass about the physical value. 

Details of simulations and the results are given in the next section, and a summary of the results in the final section. 
Details of the formalism, including expressions for various QNS are given in the appendix. 

II. SIMULATIONS AND RESULTS 
A. The simulations 

In earlier papers 0, IT^ IT^ we have shown that finite volume effects on the QNS are negligible for lattices with 
Ns > 2Nt [Ns is the spatial extent of the lattice and Nt the temporal extent) . The data we discuss here are obtained 
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FIG. 1: Some of the QNS, x/r^, as functions of T/Tc for m„d = O.lTc and m<j = T^. 



on 4 X 16'^ lattices. The setting of scale, the parameters employed and the statistics are detailed in To that 

set of data we have added two more sets — 55 configurations separated by more than two autocorrelation times at 
T/Tc = 0.975±0.010 (i.e., /3 = 5.2825) and 86 configurations, similarly spaced, at T/T^ = 1.15±0.01 (i.e., (3 = 5.325). 
The configurations are generated with a bare sea quark mass m = O.lTc, which gives = 0.3mp. 

We have explored the dependence of the physics on the strange quark mass and on variations in the light quark 
mass through partially quenched computations, i.e., the approximation in which the number of valence quark flavours 
is different from the number of dynamical sea quark flavours, and their masses are also different. Errors in partial 
quenching are bounded by comparing results with the fully quenched theory. 



B. Quark number susceptibilities 



Our primary results for QNS are shown in Figure^] These were obtained using the eqs. I|A9|I and I|A11|I in Appendix 
IXI The diagonal QNS and several of the off-diagonal ones show the characteristic crossover from small values in the 
low temperature phase to large values in the hi gh temp erature phase which gave rise to the original interpretation 
that the QCD phase transition liberates quarks [l4L Il5l|. Observe that xb < XQ through the full temperature range 
explored. Both xi and xy have values between the two others. In the low temperature phase one has xy < X/, but 
for T > 1.5Tc one obtains xy > Xi- We expect the crossover temperature between these two regimes to vary with 
quark masses. 

Our results are compatible with earlier results with staggered fermions at the same cutoff and quark mass which 
were obtained in the high temperature phase Q . They are not directly comparable to results obtained in at the 
same lattice spacing due to differences in the discretization. 



1. Robust observables 



In the quenched theory it was found that the QNS depended quadratically on the lattice spacing i.e., x(a) = 
X + 0{a'^). Since staggered Fermions have order lattice artifacts, one expects the same behaviour in the theory with 
sea quarks. We are therefore forced to search for observables which are robust against changes in the lattice spacing, 
in the sense that r{a) = r + 0{a^) with n > 2. We expect the ratios of QNS to have very good scaling properties in the 
high temperature phase, where the flavour off-diagonal QNS are much smaller than the flavour diagonal QNS. In the 
low-temperature phase we do not necessarily expect such behaviour to hold, since these two pieces are comparable, 
and the coefficient of the order corrections in the two parts depend on different physical quantities. 

As shown in Figure |21 ratios of QNS in the high temperature phase have this property. The figure also shows 
another pleasant property — these ratios have little statistically significant dependence on the sea quark content of 
the theory. We have checked that these two aspects of robustness hold for all ratios in the high temperature phase 
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FIG. 2: Ratios of QNS are robust observables — being insensitive to both changes in lattice spacing a oc 1/Nt at fixed T — 2Tc, 
and the sea quark content of QCD. The quenched results come from a reanalysis of data from 17]. In both cases the light 
valence quark mass is O.OSTc and the strange quark mass is Tc. 
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FIG. 3: Some robust predictions of fluctuation measures from QCD: all the quantities shown are the ratio Cx/q for the X 
indicated in the figure, except for X = S which is Cs/g/'i. 



of QCD. The dependence of such ratios on the valence quark masses can be determined using the quadratic response 
coefficients (QRC) defined in and apphed to the study of Cb/s- 

In view of these results, the hierarchy of QNS shown in the previous subsection must be a robust feature of QCD. 
It is therefore useful to demonstrate this hierarchy by plotting Cx/s as a function of T/Tc in Figure Our results 
indicate that experimental studies oiCg/Q, Cb/q and Cy/Q are the most promising in terms of distinguishing between 
the two phases of QCD, because they exhibit the largest changes in going from one phase to the other. 



C. Strange quarks 



The Wroblewski parameter. As, as extracted from experiments, is the ratio of the numbers of primary produced 
strange and light quark pairs. It has been argued earlier i that under certain conditions, whose satisfaction can be 
verified by independent observations, one has = Cg/^. Our results for this robust quantity are shown in Figure 01 
|19|. In this computation we have taken the strange quark mass to be rus — and the two light quark masses to 
be degenerate, m„d = O.OSTc, such that it reproduces the correct value of niT^jmp. As can be seen from the figure, 
the value of the ratio at is A^ « 0.4, in agreement with the value of the Wroblewski parameter extracted from 
experiments, when the freeze-out temperature is close to pol |. It is also a pleasant fact that at lower temperatures 
the ratio keeps decreasing. 

The dependence of this ratio on the valence quark masses was investigated in [T^. where it was shown that, in 
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FIG. 4; The robust variable C^/u = As as a function of T/Tc when the hght quark masses are taken to be rriud ~ O.OSTc, 
corresponding to a reahstic pion mass, and the strange quark mass is set to nis = Tc, which gives a reahstic value of the ratio 
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FIG. 5: The first panel shows xby /T"^ as a function of T/Tc for various patterns of SU(3) flavour symmetry breaking. Holding 
= O.lTc constant we vary nis in (A) nis — Tc, (B) — 0.75Tc, and (C) — 0.5Tc. Holding A„a = 0.25rc constant, we 
vary all the quark masses in (D) = 0.75 and (E) nis = Tc- In (F) all the quark masses are small rriud = O.OlTc, rUs = O.lTc. 
The second panel shows Xiy/T'^, as a function of T/Tc when niu ~ O.OSTc, rud — O.lTc and = Tc. 



the continuum limit, there was no dependence on valence quark mass except near T^- In the vicinity of T^,, and 
immediately below, we found Xs to be strongly dependent on mg. It increases as a function of T/Tc and at large 
enough T reaches the same value as Xui but it does this slowly when mg/T^. is large, and faster when to^ ^ T^. If 
the plasma contains strange quark quasi-particles, as we argue later, then this behaviour could be a kinematic effect, 
which measures the phase space for a thermal gluon to split into a strange quark-antiquark pair. That the first effect 
is dynamical and the second kinematical can be motivated by a study of quantities which vanish in the SU(3) flavour 
symmetric limit. 



1. Flavour symmetry breaking 

Two off-diagonal susceptibilities show an interesting pattern — Xbq and xby are both continuous through Tc, but 
peak in the vicinity of Tc. Since xbq = Xby 1'^, as seen from eqs. IjAQII and (|A11|) . we show only the latter in Figure 
for various values of quark masses explained in the caption. The figure also displays xiy for m„ 7^ m^. Direct 
computations also show that xbq = Xby — when all three quark masses are equal, and that xiy = in the SU(2) 
symmetric limit — providing an explicit demonstration that non-zero values of these quantities are due to flavour 
symmetry breaking (see the discussion in Appendix 0. 
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FIG. 6: The flavour symmetry breaking matrix elements (A) xby/^us extracted with rriud = O.OSTc and rris = O.lTc, and (B) 
^/y = X^y/^ud extracted using rriu = O.OSTc, rrid = O.lTc and = Tc, as a function of T/Tc- The kinematic suppression for 
reaUstic strange quark masses is clear from the significantly smaller values of Xsy/AL when (C) = and (D) = O.TSTc 



From a comparison of the cases (D) and (E) in FigureElit is clear that xby is not only a function of Aus — mg ~ m„ 
and T for large values of this asymmetry, since the two curves are not coincident although they have equal Aus • A 
careful look at the cases (A), (B) and (C) in the same figure shows that when nis is comparable to then both the 
position and the value of the peak in these QNS are dependent on nis- Explicit dependence of the flavour symmetry 
breaking matrix elements on the actual value of (and not just the asymmetry parameter) can only come as a 
kinematic effect. We try to confirm the magnitude of this effect next. 

In Figureiniwe display the values of the dimensionless quantities Ajy = xiy I ^'^'^ -^by = Xsy/A^,,, extracted 
using the computations in which A^s and A^d are much smaller than Tc. It would be interesting to check the 
temperature range in which these dimensionless quantities are computable in weak coupling theory. In the same 
figure we also show xsy/A^^ when A^s is comparable to T^. Its strong suppression relative to the former case shows 
the kinematic effect which is responsible for the shape of C^/u shown in Figure 01 

The physics of the region just above Tc is known to be complicated when observed through gluonic variables such 
as A/T"* — (e — 'iP)/T^ (where e is the energy density and P the pressure) as well as the ratio of the lowest lying 
screening masses in the CP-even and CP-odd sectors |23|. The peaks in Aby and Ajy are the first observations of 
interesting structures near Tc in fermionic variables unconnected with the order parameter. It would be interesting 
to see what temperature range in this is explainable by weak coupling theory. 



D. Flavour carrying degrees of freedom 



The question of which are the thermodynamically relevant degrees of freedom in the QCD plasma is easier to answer 
in the quark sector than in the gluon sector. The reason is that the multitude of flavour quantum numbers allow 
us to look for "linkage" of flavour, i.e., exciting one quantum number and seeing the magnitude of another quantum 
number that is simultaneously excited. 



1. Strangeness carriers 



Robust variables involving off-diagonal QNS serve precisely this purpose. In [2j the robust variable 

Cbs = -iC(BS)/s = -3 = 1 + — ^ = 1 + C(,us)/s + C(rfs)/s = 1 + 2C(„s)/s (4) 

is identified as one which can distinguish between bound state QCD and the usual picture of the excitations in 
the plasma phase of QCD (in the last expression above we have used eq. (IA15|I and flavour SU(2) symmetry to write 
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FIG. 7: The robust variables Cbs and Cqs, as functions of T/T^. The quark masses used are m^d ~ O.lTc and = Tc, 
although in the high temperature phase there is no statistically significant dependence on the quark masses. 



C(us)/s = C'(ds)/s)- This is expected to have a value of unity if strangeness is carried by quarks (i.e., S = 1 always 
comes linked with = — 1/3). In Q it was shown that bound state QGP gives a value of Cbs ~ 2/3 (for T > Tc). 

We present the first estimate for this quantity from lattice QCD in Figure[3 In the low-temperature phase Cbs is 
very different from unity, but immediately above Tc the value is clamped to unity. There is no statistically significant 
change in Cbs as rUs/Tc is varied between 0.1 and 1. Since the statistical error bars are extremely small for T > Tc, 
this is a strong statement which contrasts with the dependence of As and xby- 

Another interesting measure is the correlation of charge and strangeness measured by the robust observable 



When strangeness is carried by quarks one would expect this to be unity (since S — 1 comes with Q = 1/3). In Figure 
E|we have also shown the first measurement of Cqs- Immediately above Tc it reaches close to unity with small errors. 
As a result, these two measurements together quite strongly indicate that unit strangeness is carried by objects with 
baryon number —1/3 and charge +1/3 in the high temperature phase of QCD, immediately above Tc- 

Furthermore, eqs. (^ISJ indicate that our observations imply that Xus = 0, and hence strangeness carrying excita- 
tions do not carry u or d flavour. This is the most direct lattice evidence to date that strangeness is linked to other 
quantum numbers exactly as it would be for strange quarks, in the high temperature phase of QCD; and that these 
linkages are quite different below Tc- Later in this section we show that one should think of these as quasi-particles, 
dressed by the strong residual interactions, rather than as elementary quarks. 

Apart from the direct evidence of linkage between quantum numbers, we also draw attention to the cryptic evidence 
in the temperature and dependence of Cbs and Cqs- The rapid change of Cbs with T (for T < Tc) has a natural 
explanation if the thermodynamics is controlled by a spectrum of strange baryons such that the amount of (anti-) 
strangeness per baryon increases with mass, and the masses are larger than T. The temperature independence of 
the two quantities above Tc similarly implies that there is one excitation, which has mass less than Tc- The fact 
that the values of these quantities do not depend on within errors, for T > Tc, further implies that the effective 
masses of these quasi-particles is less than Tc, so that the infrared cutoff on the Dirac operator spectrum is provided 
by T. The next heavier quark, charm quark, does not affect the thermodynamics of the QCD plasma in the range 
of temperatures we investigate, since its mass is well beyond T. However, these heavier quarks do probe changes in 
other aspects of physics, such as screening, as is evident from [23| . 



Cqs - 3C(Qs)/5 - 1 



2Xtis — Xds 



(5) 



2. The light quark sector 



In transplanting these methods to the light quark sector, we find that the composite QNS, xbi and XQi^ are not 
informative, since the quark of one flavour has the same isospin as the antiquark of the other flavour. One way to 
extract information on the degrees of freedom would be to consider QNS of G-parity. However, it is more transparent 
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FIG. 8: The robust variable ~C(^ud)/u which measures the correlation between u and d flavours for rriud = O-lTc. It is positive 
in the low temperature phase since u quarks are found along with d antiquarks in charged pions, and vanishes in the high 
temperature phase, indicating that u and d are fully decorrelated in the plasma. There is no statistically significant dependence 
on niud in the high temeperature phase. 



to turn to the flavoured QNS Xud oc {MuMd)- We can then use the quantity 

C(ud)/u = (6) 

Xu 

which looks at the hnkage between u and d flavours in the same way that C(^qs)/s looked for linkage of strangeness 
and charge. Our results are plotted in the first panel of Figure |S1 In the hadronic phase it is non- vanishing because 
of charged pions, and negative because in these mesons each u comes with a d, and vice versa. In the QGP phase the 
vanishing of this normalized covariance implies that a particle with u quantum number does not exhibit d quantum 
numbers. 

Further tests come from investigating 

C(BU)/U = C{BD)/D = {ud)/u ~^ C!{us) /u) J 

C(QU)/U = 2 (2 — C(ud)/.u — C(i,s)/„), 

C{QD)/D = - 2C(„d)/ti + ^(us)/^)- (7) 

The vanishingly small values of C(^ud)/u ^nd C(„s)/u imply that the u flavour is carried by excitations with baryon 
number +1/3 and charge +2/3, whereas the d flavour is carried by particles with baryon number +1/3 and charge 
— 1/3. These are, therefore, quark quasi-particles. 



3. Quasi-quarks 



One might wonder why we talk of, for example, baryon number 1/3 when the measurements even at 2Tc differ 
from this number by a few parts in a thousand. What does this small but statistically significant deviation tell 
us? The answer is that it says something about the spatial structure of the quasi-particle. If flavour were carried 
by pointlike bare quarks, then Xud and Xus would be precisely zero. However, interactions dress each quark into a 
spatially extended quasi-particle, and a thermodynamic average probes the spatial dimension of the charge with a 
resolution of 1/2-kT. When T is sufficiently large, so that the gauge coupling is sufficiently small, this structure can 
be computed in weak coupling theory. As the coupling grows, the perturbative computation fails quantitatively, but 
as long as the correction to the charge or baryon number remains small, one can fruitfully talk of quasi-quarks. 

In Figure|n|we show the flavour off-diagonal QNS Xud/T"^ for two different quark masses, along with the prediction 
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FIG. 9: The off-diagonal QNS Xud/T'^ for two different quark masses compared to weak coupling theory. The band includes 
uncertainties due to neglected of higher loop effects, the effect of changing from one-loop to two- loop computation of the running 
coupling, and statistical uncertainties in the determination of Tc/A-jjg. 



of weak coupling perturbation theory — 

Xud 10 3 1 c \ 

where c is a constant whose evaluation requires larger number of loops in the perturbation theory. The strong coupling, 
as, has been evaluated to two loop accuracy at scale 2ttT with the estimate Tc/Ajj^ = 0.49 ±0.05 0. This variation 
in Tc/A-jTs, a variation of c by two orders of magnitude, 0.1 < c < 10, and the variation in as in going from one loop to 
the two loop expression are included in the band in the figure. We find that in this range of temperature the prediction 
is somewhat smaller than the lattice data. Since this is not a robust variable, it is possible that taking the continuum 
limit will improve the agreement between the two. However, it is clear that as one comes closer to Tc the disagreement 
increases, although the magnitude of C(„d)/ti remains small. Thus, it seems that a Fermi gas picture, which may be 
valid at large T/Tc, gives way to something more complicated as one approaches Tc, although the quantum numbers 
are linked in exactly the same way as for the elementary quarks. This is the meaning of quasi-quarks. 



III. SUMMARY 



We have presented an extensive computation of many different quark number susceptibilities (see Appendix ^ for 
the definitions). All the diagonal QNS, and some of the off-diagonal QNS, track the phase structure of QCD — being 
small in the confined phase and crossing over to larger values in the high temperature phase of QCD, as shown in 
Figure ^ 

An important observation was that ratios of QNS, Ca/b, defined in eq. are robust variables which depend 
weakly on the lattice spacing and the sea quark content of QCD in the high temperature phase, as shown in Figure 
121 These ratios can be compared to experimentally determined ratios of variances (or covariances) in event-to-event 
fluctuations of conserved quantum numbers. The relative magnitudes of the diagonal QNS are among these robust 
observables, and we found the ordering xs > Xq > Xy > Xi > Xb, shown in Figure |21 

A second set of results concerns the thermal production rate of strange quarks. It has been argued [T^ that under 
certain (testable) conditions the Wroblewski parameter is the robust observable Cs/u- While it is insensitive to the sea 
quark content of QCD, it is known to depend sensitively on the valence quark masses [T^ . Here we have determined 
this quantity for realistic values of the strange and light quark masses (see Figure 0}. 

We attributed this dependence on nis to kinematic effects visible when nis ~ 0{Tc). However, kinematic effects 
should manifest themselves in other quantities as well. We tested this hypothesis by examining certain QNS which 
vanish in the flavour symmetric limit. We extracted the matrix elements which are quadratic in the flavour symmetry 
breaking mass differences, A„s, when A„s <gC T^.. By comparing these (in Figure I^Jl to the corresponding quantities 
when Aus ~ OiT^), we demonstrated the presence of such kinematic effects in other quantities as well. The flavour 
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symmetry breaking matrix elements themselves (Figure peak at T slighly larger than Tc, and are the first known 
example of observables in the quark sector of QCD which parallel similiar structures seen in the gluon sector. 

Our final result is that the high temperature phase of QCD essentially consists of quasi-quarks. We demonstrated 
this by observing that unit strangeness is carried by something which has baryon number —1/3 and charge 1/3, as 
shown in Figure[7| Part of the argument is that this correlation does not depend on the strange quark mass even when 
it is as large as Tc. Similarly, in the light quark sector one finds that u and d quantum numbers are not produced 
together (Figure (S)). Through eqs. Q we found that this implies that the u flavour is carried by excitations with 
baryon number +1/3 and charge +2/3, whereas the d flavour is carried by particles with baryon number +1/3 and 
charge —1/3. 

We presented an argument that the carriers of these quantum numbers are not elementary quarks but their dressed 
counterparts which are called quasi-quarks. This argument involved the comparison of Xud with a weak-coupling 
prediction, which is shown in Figure |^ The key point is that this comparison fails badly as one approaches Tc, 
although the correlations of flavour quantum numbers remains as they would for quarks. A similar comparison of 
weak coupling prediction with lattice results for the diagonal QNS Xu also fails near Tc, leading us to the same 
conclusion. 

The argument about the existence of quasi-quarks in the high temperature phase of QCD depends on the exami- 
nation of robust variables given in eq. Q. It is useful to note that their use is not restricted to the lattice. It is also 
possible to measure them in heavy-ion collisions and thereby deduce the nature of excitations in the fireball produced 
in these collisions. 

We end by pointing out that we have not studied the low-temperature phase of QCD in much detail here. This is 
an interesting problem, which we have touched upon very briefly in the discussion oi Cbs and Cqs, and has been left 
for the future. 
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on the Cray. Part of this work was done during a visit under Indo-French (IFCPAR) project, 3104-3 to SPhT, Saclay. 
The hospitality of Saclay and IFCPAR's support is gratefully acknowledged. 



APPENDIX A: SUSCEPTIBILITIES IN EQUIVALENT ENSEMBLES 

Corresponding to every conserved charge, Q, under the global symmetries of a theory, one can introduce a chemical 
potential, /i, into thermodynamics, by adding to the action a source term ^Q. In QCD, at finite quark mass, one has 
SU{Nf) vector flavour symmetry. Corresponding to each of the Nj flavours, one can introduce a chemical potential 
fif (/ = u, d, s, etc.) through the term 

J = Y,t,fMf = ^i^M, (Al) 
/ 

where J\ff is the number operator for quarks of flavour /, whose expectation value is the number of quarks minus the 
number of antiquarks. The last expression just rewrites the sum as a dot product of the vector of intensive variables 
fi with the vector of extensive variables M. This corresponds to a grand canonical ensemble in which the chemical 
potential on each of the quark flavours can be tuned independently. In the corresponding canonical ensemble the 
number of quarks of each flavour is kept fixed. 

The number densities, which are the first derivative of the pressure with respect to the chemical potential, and the 
quark number susceptibilities (QNS), which are the second derivatives, have been defined before (1^, 17J. Here we use 
the notation of 0| for the QNS. We shall also use a higher order susceptibility, for which we use the notation of [T^ . 

It is usually more convenient to define chemical potentials for variables which are easier to control in experiments — 
such as the baryon number, B, the electric charge, Q, the third component of isospin, I3, or the hypercharge, Y. Any 
choice of Nf variables corresponds to a different choice of ensemble to describe the same physics. The description in 
terms of fiavours given above can then be translated into the new ensemble by a simple linear transformation 

J = n^M-HlAf = Ui'fAf', i.e., TV' = MJ\f and ^i' = {M-^f^i = {M^y^^i. (A2) 

Clearly, the partition function being the same, the physics remains invariant under these redefinitions. The choice of 
a given M corresponds to putting coordinates in Gibbs space. 

One is usually interested in thermodynamics quantities or response functions which are obtained by taking deriva- 
tives of the free energy or pressure with respect to the chemical potentials. Note that by the definitions in eq. I|A2|I . 
one has /i = M'^ ^' . The chain rule for differentiation then tells us that 

A_^A-fM^V — -M— (A3) 
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The fact that this gives back the original definitions of the Af's in terms of the Afs is a consistency check of the 
formahsm. We illustrate the uses of this formalism for the cases of Nf = 2 and Nj = 3 below. 

1. Nf = 2 

a. The B, I3 ensemble 

In the two flavour case, one can transform from the flavour basis to the set 

1 3 

h = ^(A/; - Nd), m = (A4) 

Inverting the relation between chemical potentials one obtains /i„ — (Ib/S + fJ-i/'2 and fid — (J-b/'^ — In the 

ensemble where nj = 0, one then gets = fJ-d — fJ'B/'i- The number densities are given by eq. (|A4() . The quark 
number susceptibilities are — 

1 2 

XB = g(x«« +Xd<i + 2x„d) = g(x« +Xud), 

XI = -^{Xuu + Xdd - 2Xud) = -j^iXu - Xitd) , 

XBi = ^{Xuu - Xdd) = 0, (A5) 

where the first expression in each case is the most general, and the second is obtained for exact vector SU(2) flavour 
symmetry m^d- If this symmetry is broken then xbi should become non-zero. For small values of the symmetry 
breaking parameter A^d = — m„ (we will follow the convention that m„ < < ms), 

XBi^ABi^la, which yields ^ ^ Abi 1^^^ , (A6) 

where Agj \s a, dimensionless number. We present results for this quantity in Figure |S1 In the low temperature phase 
we expect that Abi is a non-perturbative quantity, but that it should be computable in chiral perturbation theory. 
It would be interesting to check how far the weak coupling theory in the high temperature phase agrees with our 
determination of Abi- 

b. The B, Q ensemble 
One may choose to work in another ensemble given by 

B = ^(A^u + Afd), fiB = + 2^/^, 

Q = ^(2A^„ - Afd), fJ-Q = fJ-u - fJ-d- 

For HQ = 0, one gets again the expected result = i^d = I^b/^- The quark number susceptibilities in this ensemble 
are — 

1 2 

Xb = g (x«« + Xdd + 2xu(i) = g(x« + x«d)j 

XQ = ^(4x™ + Xdd-4x„d) = ^(5xn-4x«d), 

XBQ = ^(2x™-Xdd + Xud) = ^(X« + X«d), (A7) 

where the last expression in each line holds in the special case of ruud- Note that xbb is the same in both the 
ensembles. This follows from the fact that the definition of the baryon number is the same. 
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2. iV^ = 3 

a. The B, I^, Y ensemble 



The variables in this ensemble are 



-^3 = ^(A^u - ^fd), ^^I = M« - Mrf, 

^(A/'u+A/'<i-2A/;) My = ^(M« + Md-2Ms)- (A8) 



The six independent quark number susceptibilities are 



XB = 


g iXu 


+ Xd + 


Xs + 2x«d + 2x«s + 2x<is) = 


= ^(2X« 


+ Xs + 2x«d + 


4x«s), 


XI = 


■^{Xu 


+ Xd- 


2XtKi) = 2^Xu- Xud), 








XY = 


g iXu 


+ Xd + 


4Xs + 2x«d - 4:Xus ~ 4xds) 


~ g(XM 


+ 2xs + Xud - 


4Xtis)j 


XBI = 


g iXu 


-Xd + 


Xus - Xds) = 0, 








Xby = 


g (X« 


+ Xd~ 


2Xs + 2x«d - Xus - Xds) ^ 


^(Xn- 


Xs Xud Xui 




XIY = 


g (X« 


-Xd~ 


2xus + 2xds) = 0. 









(A9) 

As before, the last set of expressions on each line holds only for niud 7^ ^s- Similiar to eq. (|A6I) . one can define Abi 
here, and also Ajy — xiy / ^td^ both of which are generally non- vanishing when SU(2) flavour symmetry is broken. 

In the SU(3) symmetric limit, m„ — — rus, the three off-diagonal susceptibilities vanish, i.e., xbi = Xby = 
XiY — 0. Also, in this limit one has 2x/ — 3xf/2 ~ Xu — Xud and 3xb — Xu + '^Xud- In the low temperature phase the 
breaking of vector SU(3) symmetry produces the mass difference between the pion and the K meson. At sufficiently 
high temperature, when the strange quark mass is less than the Matsubara frequency, < 27rT, the theory becomes 
effectively SU(3) symmetric, and the above relations should hold. In the high temperature phase of QCD one also 
has Xud g^logg ^0 @, so one should obtain 2x/ = 3xy/2 = 3xb- 

b. The B, Q, Y ensemble 
Another useful set of charges and associated chemical potentials is 

B = i(7V„ +Nd+ Ms), Ms = Mu + Md + Ms, 

Q = i (27V„ -JVd-Ms), ^ flu- fJ'd, 

Y = r^{Mu+Ki-2J\fs), fiY^fid-fis- (AlO) 

The three susceptibilities xs, Xy and Xby f^i'e the same as in the previous ensemble. The remaining susceptibilities 
are 

XQ = ^(^Xu + Xd + Xs- ^Xud - 4Xus + 2xds) = ^(5x« + Xs - 4xud - 2x„s), 
XSQ = n(2x " 2Xds) = g(Xu - Xs + Xud - X«s), 

Xqy = ^(2Xu-Xd + 2xs+X«d-5x«s+Xds) = ^(Xu + 2xs + Xud-4x«s)- (All) 

As before, the last set of expressions on each line holds only for m„(i ^ mg. Note that in this limit one has xbq — 
XbyI2. 
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In the SU(3) symmetric limit, m„d = m^, two of the off-diagonal susceptibilities vanish, i.e., xbq = Xby = 0. Also, 
in this hmit one has xq ^ Xy = "^XQy = 2(x« - Xud)/^ and xb = {Xu + 2x«<i)/3. As before, in the high temperature 
limit, when the quark masses are less than the temperature, the theory becomes effectively SU(3) symmetric. Then 
taking Xud = 0, one should obtain xq = Xy = '^Xb = ^XQY- 

It is interesting to examine this in a theory of massless free fermions. The free energy is given by 



F = NeV J2 

f—u.d,s 



180 +6^^ ^ 12n^^ 



(A12) 



Substituting the values of the flavour chemical potential by the appropriate combination of mb, Mq and /xy, and 
taking the derivatives, we find that xby = Xbq = 0. Also, XQ = Xy = "^Xb ~ "^Xqy = 2Nc/9. For massive free 
fermions, when T is much larger than the fermion mass, the same results would hold. 

When SU(3) symmetry is broken through the parameter A„s = rus — to„ (we assume that SU(2) symmetry still 
holds) then for small A^s one may again write 

Xby = AbyAI^ and xbq = AbqAI^^- (A13) 

As before, we expect Aby and Abq to be non-perturbative but computable in chiral perturbation theory in the low 
temperature phase. It would be interesting to compare our results (presented later) with weak coupling theory in the 
high temperature phase. 

c. The B, Q, S ensemble 
From the experimental point of view, it may be interesting to use the set 

B = ^i^fu+^fd+^fs), ij,b = tJ.u + '^iJ'd, 

Q = ^(2A/'„ -Nd- Ms), mq = Mtt - Md, 

S = -Ns, ns = tM-t^s- (A14) 

Note that we have used the standard convention where the strangeness of the antistrange quark is +1. The three 
susceptibilities, xb, XQ and xbq are as before. The remainder are 

Xs = Xs, 

XBS = -^{Xa + Xua+Xds) = -^{Xs + '^Xus)-, 

Xqs = ^{Xs - 2Xus + Xds) = ^(Xs - Xus)- (A15) 
As always, the last set of expressions on each line holds only for m„d ^ m^. 

d. The B, Q, U ensemble 

For technical questions about the light quark sector it is useful to work in the ensemble with 

Q = l{2Nu-Md-K), 

U = Mu, (A16) 

The three susceptibilities, xb, XQ and xbq are as before. The rest are 

Xu = Xu, 
1 

Xbu — 2 (X« "I" Xud "I" Xus ) ) 

XQU = \{'^Xu-Xud-Xus)- (A17) 
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Changing to an ensemble where U is replaced by D = Md changes the QNS to 

XD = Xu, 
Xbd = -^{Xu + Xud + Xus) , 

XQD = -^{Xu-2xud + Xus), (A18) 
where we have used SU(2) symmetry. This gives xu = Xd and xbu — Xbd- 
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